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Coronaviruses are enveloped, positive-stranded RNA viruses with a genome of approximately 30 kb. Based 
on genetic similarities, coronaviruses are classified into three groups. Two group 2 coronaviruses, human 
coronavirus OC43 (HCoV-OC43) and bovine coronavirus (BCoV), show remarkable antigenic and genetic 
similarities. In this study, we report the first complete genome sequence (30,738 nucleotides) of the prototype 
HCoV-OC43 strain (ATCC VR759). Complete genome and open reading frame (ORF) analyses were performed 
in comparison to the BCoV genome. In the region between the spike and membrane protein genes, a 290- 
nucleotide deletion is present, corresponding to the absence of BCoV ORFs ns4.9 and ns4.8. Nucleotide and 
amino acid similarity percentages were determined for the major HCoV-OC43 ORFs and for those of other 
group 2 coronaviruses. The highest degree of similarity is demonstrated between HCoV-OC43 and BCoV in all 
ORFs with the exception of the E gene. Molecular clock analysis of the spike gene sequences of BCoV and 
HCoV-OC43 suggests a relatively recent zoonotic transmission event and dates their most recent common 
ancestor to around 1890. An evolutionary rate in the order of 4 x 10 -4 nucleotide changes per site per year was 
estimated. This is the first animal-human zoonotic pair of coronaviruses that can be analyzed in order to gain 
insights into the processes of adaptation of a nonhuman coronavirus to a human host, which is important for 
understanding the interspecies transmission events that led to the origin of the severe acute respiratory 
syndrome outbreak. 


Coronaviruses are large (120- to 160-nm), roughly spherical 
particles with a linear, nonsegmented, capped, and polyade- 
nylated positive-sense single-stranded RNA genome that is 
encapsidated in a helical nucleocapsid. The envelope is derived 
from intracellular membranes and contains a characteristic 
crown of widely spaced club-shaped spikes that are 12 to 24 nm 
long. The genus Coronavirus (International Committee on the 
Taxonomy of Viruses database [ICTVdb], virus code 
03.019.0.1) belongs to the family Coronaviridae in the order 
Nidovirales (7, 8). 

Before the 2002-to-2003 severe acute respiratory syndrome 
(SARS) epidemic, coronaviruses were somewhat neglected in 
human medicine, but they have always been of considerable 
importance in animal health. Coronaviruses infect a variety of 
livestock, poultry, and companion animals, in whom they can 
cause serious and often fatal respiratory, enteric, cardiovascu¬ 
lar, and neurologic diseases (25). Most of our understanding 
about the molecular pathogenic properties of coronaviruses 
has been achieved by the veterinary virology community. 

The coronaviruses are classified into three groups based on 
genetic and serological relationships (19). Group 1 contains 
the porcine epidemic diarrhea virus (PEDV), porcine trans¬ 
missible gastroenteritis virus (TGEV), canine coronavirus 
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(CCoV), feline infectious peritonitis virus (FIPV), human 
coronavirus 229E (HCoV-229E), and the recently identified 
human coronavirus NL63 (HCoV-NL63). Group 2 contains 
the murine hepatitis virus (MHV), bovine coronavirus 
(BCoV), human coronavirus OC43 (HCoV-OC43), rat sialo- 
dacryoadenitis virus (SDAV), porcine hemagglutinating en¬ 
cephalomyelitis virus (PHEV), canine respiratory coronavirus 
(CRCoV), and equine coronavirus (ECoV). Group 3 contains 
the avian infectious bronchitis virus (IBV) and turkey corona¬ 
virus (TCoV). The SARS coronavirus (SARS-CoV) is not as¬ 
signed to any of these groups but is most closely related to 
group 2 coronaviruses (21, 54). 

HCoV-OC43 (ICTVdb code 19.0.1.0.006) and HCoV-229E 
(ICTVdb code 19.0.1.0.005) were isolated in 1967 from volun¬ 
teers at the Common Cold Unit in Salisbury, United Kingdom. 
HCoV-OC43 was initially propagated on ciliated human em¬ 
bryonic tracheal and nasal organ cultures (42). HCoV-OC43 
and HCoV-229E are responsible for 10 to 30% of all common 
colds, and infections occur mainly during the winter and early 
spring (38). The incubation period is 2 to 4 days. During the 
2002-to-2003 winter season, a new human coronavirus, HCoV- 
NL63, was isolated from a 7-month-old child suffering from 
bronchiolitis and conjunctivitis in The Netherlands (61). Seven 
additional HCoV-NL63-infected individuals, both infants and 
adults, were identified, indicating that HCoV-NL63 can be 
considered an important new etiologic agent in respiratory 
tract infections. Coronaviruses infect all age groups, and rein¬ 
fections are common. The infection can be subclinical and is 
usually mild, but there have been reports of more-severe lower 
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respiratory tract involvement in infants and elderly people (17, 
60). Human coronaviruses can induce a demyelinating disease 
in rodents and can infect primary cultures of human astrocytes 
and microglia. A possible etiological role for HCoV-OC43 and 
HCoV-229E in multiple sclerosis is being debated (4, 13, 15). 

The coronavirus genomes are the largest of the known RNA 
viruses (27 to 31.5 kb) and are polycistronic, generating a 
nested set of subgenomic RNAs with common 5' and 3' se¬ 
quences (35). The 5' two-thirds of the genome consists of two 
large replicase open reading frames (ORFs), ORFla and 
ORFlb. The ORFla polyprotein (ppla) can be extended with 
ORFlb-encoded sequences via a —1 ribosomal frameshift at a 
conserved slippery site (6), generating the >7,000-amino-acid 
polyprotein pplab, which includes the putative RNA-depen- 
dent RNA polymerase (RdRp) and RNA helicase (HEL) ac¬ 
tivity (20, 39). The polyproteins ppla and pplab are autocata- 
lytically processed by two or three different viral proteases 
encoded by ORFla: one or two papain-like proteases (PLP1 
and PLP2) and a 3C-like protease (3CL pro ) (39, 67, 68). Other 
putative domains presumably associated with a 3'-to-5' exonu¬ 
clease (ExoN) activity, a poly(U)-specific endo-RNase (Xen- 
doU) activity, and a 2'-0-methyltransferase (2'-0-MT) activity 
are predicted in pplab (27, 54). The 3' end of a coronavirus 
genome includes several structural and accessory protein 
genes: an envelope-associated hemagglutinin esterase (HE) 
glycoprotein gene, present only in group 2 coronaviruses; a 
spike (S) glycoprotein gene; an envelope (E) protein gene; a 
matrix (M) glycoprotein gene; a nucleocapsid (N) phospho- 
protein gene; and several ORFs that encode putative nonstruc- 
tural (ns) proteins (35). 

Coronaviruses are well equipped to adapt rapidly to chang¬ 
ing ecological niches by the high mutation rate of their RNA 
genome (about 10~ 4 nucleotide substitution/site/year) and 
high recombination frequencies (51). Many animal coronavi¬ 
ruses cause long-term or persistent enzootic infections. Long 
periods of coronavirus infection combined with a high muta¬ 
tion and recombination rate increase the probability that a 
virus mutant with an extended host range might arise. 

The current emergence of the SARS-CoV is an example of 
a crossing of the animal-human species barrier. It is likely that 
the SARS-CoV was enzootic in an unknown animal or bird 
species before suddenly emerging as a virulent virus for hu¬ 
mans. Chinese scientists found that six masked palm civets 
(Pagiima larvata ) and a racoon dog (Nyctereutes procyonoides) 
for sale in an exotic food market in Shenzhen, in the Guang¬ 
dong province in Southern China, were harboring a virus very 
similar to the SARS-CoV (1). Thirteen percent of the civet 
merchants tested at markets in Guangdong also had SARS 
antibodies. Sequence analysis showed that the animal version 
of the SARS-CoV contained an extra stretch of 29 bases (22). 
It is still not clear whether the civets were a reservoir for the 
virus or were infected by another species. 

HCoV-OC43 and BCoV (ICTVdb code 03.019.0.01.002) 
show remarkable antigenic and genetic similarities (23, 29, 36, 
44, 52, 63, 65). They both have hemagglutinating activity by 
attaching to the 7V-acetyl-9-0-acetylneuraminic acid moiety on 
red blood cells (33). BCoV causes severe diarrhea in newborn 
calves. The complete nucleotide sequences of different BCoV 
strains are known, but only fragments of the HCoV-OC43 
genome had been determined previously. In this paper, we 


report the complete HCoV-OC43 sequence (30,738 bases) and 
the comparative characterization and evolutionary relationship 
of the BCoV-HCoV-OC43 pair. This is the first animal-human 
zoonotic pair of coronaviruses that can be analyzed in order to 
gain insights into the processes of adaptation of a nonhuman 
coronavirus to a human host. 

MATERIALS AND METHODS 

Preparation of HCoV-OC43 RNA. An HCoV-OC43 strain (VR759) was ob¬ 
tained from the American Type Culture Collection (ATCC). The ATCC VR759 
strain originated from a volunteer with a common cold-like illness at the Com¬ 
mon Cold Unit in Salisbury, United Kingdom (42). HCoV-OC43 was propagated 
in a human rhabdomyosarcoma (RD) cell line, obtained from the European 
Collection of Cell Cultures (ECACC). The supernatant was harvested after 7 
days of incubation at 33°C, and RNA was isolated by using the QIAamp viral 
RNA kit (QIAGEN, Westburg, The Netherlands). A real-time quantitative re¬ 
verse transcription PCR (RT-PCR) (Taqman; Perkin-Elmer Applied Biosys¬ 
tems, Foster City, Calif.) was developed to determine the number of RNA copies 
present in the supernatant. 

Sequencing of the HCoV-OC43 genome. To determine the HCoV-OC43 
genomic sequence, a set of overlapping RT-PCR products (average size, 1.5 kb) 
encompassing the entire genome was generated. For both RT-PCR and sequenc¬ 
ing, oligonucleotide primers were designed in regions that were conserved be¬ 
tween the BCoV and MHV genomes. The forward PCR primer in the 5'- 
terminal sequence (OC43F1 [5'-GATTGTGAGCGATTTGC-3']) was based on 
the HCoV-OC43 5' untranslated region partial sequence (H. Y. Wu, J. S. Guy, 
D. Yoo, R. Vlasak, and D. A. Brian, unpublished data; GenBank accession 
number AF523847). To generate RT-PCR products containing the exact 3'- 
terminal sequence, we used oligonucleotide OC43R74 (5'-TTTTTTTTTTGTG 
ATTCTTCCA-3') based on the conserved 3'-end sequence of all known group 
2 coronaviruses. By using 150 sequencing primers, sequencing in both directions 
was performed on an ABI Prism 3100 genetic analyzer (Perkin-Elmer Applied 
Biosystems) using the BigDye terminator cycle sequencing kit (version 3.1). 
Chromatogram sequencing files were inspected with Chromas 2.2 (Technely- 
sium, Helensvale, Australia), and contigs were prepared by using SeqMan II 
(DNASTAR, Madison, Wis.). 

DNA and protein sequence analyses. ORF analysis was performed by using the 
NCBI ORF finder (http://www.ncbi.nlm.nih.gov/gorf/gorf.html). Potential SC- 
like protease cleavage sites were identified by using the NetCorona 1.0 server 
(30). DNA and protein similarity searches were performed using the NCBI 
WWW-BLAST (basic local alignment search tool) server on the GenBank DNA 
database, release 118.0 (2). Pairwise nucleotide and protein sequence alignments 
were performed by using FASTA algorithms in the ALIGN program on the 
GENESTREAM network server (http://www2.igh.cnrs.fr) at the Institut de Ge- 
netique Humaine in Montpellier, France (47). Maizel-Lenk dot matrix plots 
were calculated using the pairwise FLAG 1.0 (fast local alignment for gigabases) 
algorithm at the server of the Biomedical Engineering Center of the Industrial 
Technology Research Institute in Hsinchu City, Taiwan (http://bioinformatics 
.itri.org.tw/prflag/prflag.php). Multiple sequence alignments were prepared by 
using CLUSTALW (58) and CLUSTALX, version 1.82 (59) and were manually 
edited in GENEDOC (46). Phylogenetic analyses were conducted by using 
MEGA, version 2.1 (34). 

Evolutionary rate analyses and timing of the most recent common ancestor. 

The relationship between isolation date and genetic divergence was investigated 
using a linear regression, based on a maximum-likelihood tree, as implemented 
in the Path-O-Gen software, kindly provided by Andrew Rambaut (University of 
Oxford, Oxford, United Kingdom). Evolutionary rates and divergence times 
were estimated by using maximum likelihood in the TipDate software package, 
version 1.2 (49), and Bayesian inference in BEAST, version 1.03 (kindly made 
available by A. J. Drummond and A. Rambaut, University of Oxford; http: 
//evolve.zoo.ox.ac.uk/beast/). The molecular clock hypothesis was tested by using 
the likelihood ratio test. 

Nucleotide sequence accession number. The nucleotide sequence data reported 
in this paper were deposited in GenBank under accession number AY391777 by 
using the National Center for Biotechnology Information (NCBI; Bethesda, 
Md.) Banklt v3.0 submission tool (http://www3.ncbi.nlm.nih.gov/BankIt/). 

RESULTS 

HCoV-OC43 complete genomic sequence. We report here 
the complete nucleotide sequence of the prototype HCoV- 
OC43 strain (VR759), isolated in 1967 from an adult with 
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FIG. 1. Linear representation of the ORFs of the group 2 coronaviruses and SARS-CoV. Nucleotide insertions (open arrowheads) and 
deletions (solid arrowheads) in the HCoV-OC43 genome compared to BCoV are shown. 


common cold-like symptoms (42). The HCoV-OC43 genome 
encompasses 30,738 nucleotides [excluding the 3' poly(A) tail] 
and was deposited in the GenBank database under accession 
number AY391777. The HCoV-OC43 genome has a GC-con- 
tent of 36.9%. 


ORF organization of HCoV-OC43. The HCoV-OC43 ge¬ 
nome contains 11 major ORFs flanked by 5' and 3' untrans¬ 
lated regions of 211 and 288 nucleotides, respectively. A linear 
representation of the major ORFs of FICoV-OC43, other 
group 2 coronaviruses, and SARS-CoV is given in Fig. 1. Table 
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TABLE 1. Positions of the major ORFs of HCoV-OC43 (ATCC VR759) and BCoV (Mebus strain) 


ORF 


Nucleotide position 


No. of 
bases 

No. of 
amino acids 

Virus 

Start ORF 

First ATG 

Stop codon 

ORFla 

HCoV-OC43 

205 

211 

13362 

13,158 

4,385 


BCoV 

205 

211 

13362 

13,158 

4,385 

ORFlb 

HCoV-OC43 

13332 

13566 

21497 

8,166 

2,721 


BCoV 

13332 

13566 

21494 

8,163 

2,720 

ns2 

HCoV-OC43 

21498 

21507 

22343 

846 

281 


BCoV 

21495 

21504 

22340 

846 

281 

HE 

HCoV-OC43 

22334 

22355 

23629 

1,296 

431 


BCoV 

22331 

22352 

23626 

1,296 

431 

S 

HCoV-OC43 

23623 

23644 

27729 

4,107 

1,368 


BCoV 

23620 

23641 

27732 

4,113 

1,370 

ns4.9 

HCoV-OC43 

NA“ 

NA 

NA 

NA 

NA 


BCoV 

27811 

27889 

28026 

216 

71 

nsl2.9 

HCoV-OC43 

27802 

27817 

28146 

345 

114 


BCoV 

28095 

28110 

28439 

345 

114 

E 

HCoV-OC43 

28133 

28133 

28387 

255 

84 


BCoV 

28426 

28426 

28680 

255 

84 

M 

HCoV-OC43 

28381 

28402 

29094 

714 

237 


BCoV 

28674 

28695 

29387 

714 

237 

N 

HCoV-OC43 

29095 

29104 

30450 

1,356 

451 


BCoV 

29388 

29397 

30743 

1,356 

451 

la 

HCoV-OC43 

29102 

29165 

29347 

246 

81 

lb 

HCoV-OC43 

29348 

29441 

29788 

441 

146 

I 

BCoV 

29395 

29458 

30081 

687 

228 


a NA, not applicable. 


1 shows a comparison of the positions of the major ORFs of 
HCoV-OC43 and BCoV strain Mebus. 

The first two-thirds of the genome consists of two large 
replicase ORFs, ORFla and ORFlb. ORFla is 4,383 codons 
long and overlaps with ORFlb, which consists of 2,721 codons. 
The coronavirus replicase polyprotein of 7,095 amino acids 
(aa) is synthesized by a — 1 ribosomal frameshift at a conserved 
slippery site (UUUAAAC, nucleotides 13335 to 13341). In this 
polyprotein (GenBank accession number AAR01012) numer¬ 
ous putative functional domains are predicted: PLP1 and PLP2 
(aa 852 to 2750), 3CL pro (aa 3247 to 3549), RdRp (aa 4370 to 
5297), HEL (aa 5298 to 5900) (20, 39), and several putative 
nidovirus homologs of cellular RNA-processing enzymes, such 


pplab residues: 

V V T ▼ 


PLP1, X, PLP2, TM1 || TM2 

^RnnnR 

MB,HEL, 
ATPase 

exoN 

XendoU 

- 


i-1-1-1-1-1-1-r 

1 1000 1000 1000 0000 0000 6000 7000 

FIG. 2. Overview of the putative domain organization and potential proteolytic cleavage sites of the HCoV-OC43 replicase polyprotein pplab. 
Cleavage sites that are predicted to be processed by 3C-like protease are indicated by black arrowheads, while potential papain-like protease 
cleavage sites are indicated by white arrowheads. The following predicted domains are shown: papain-like proteases 1 and 2 (PLP1 and PLP2), 
X domain (X), putative transmembrane domains 1, 2, and 3 (TM1, TM2, and TM3), 3C-like protease (3CL), growth factor-like domain (GFL), 
RdRp, metal ion-binding domain (MB), HEL, ATPase, putative 3'-to-5' exonuclease (ExoN), putative poly(U)-specific endo-RNase (XendoU), 
and a putative 5-adenosylmethionine-dependent ribose 2'-0-methyltransferase (MT). 


as ExoN (aa 5901 to 6421), XendoU (aa 6422 to 6796), and a 
2'-0-MT domain (aa 6797 to 7095) (Fig. 2) (27, 54). The 
amino-terminal part of the polyprotein is predicted to be 
cleaved by the papain-like proteases (68), while the carboxy- 
terminal part is putatively processed by the main coronavirus 
protease, 3CLP ro (67). 

The 3'-proximal part of the HCoV-OC43 genome contains 
several ORFs, which encode a variety of structural and acces¬ 
sory proteins. Downstream of ORFlb, a nonstructural protein 
gene (ns2) of 837 nucleotides, which is a group 2-specific gene, 
is present. Although these group 2-specific genes are not es¬ 
sential for viral growth, recent work has shown that deletion of 
MHV ns2 leads to a significant attenuation of the virus when 
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TABLE 2. Nucleotide and amino acid similarities of the major HCoV-OC43 (ATCC VR759) ORFs with the ORFs of BCoV, CRCoV, 

PHEV, ECoV, MHV, and SDAV 


HCoV-OC43 

ORF 



% Nucleotide (amino acid) similarity 



BCoV 

CRCoV 

PHEV 

ECoV 

MHV 

SDAV 

ORFla 

97.4 (97.0) 

NA" 

NA 

NA 

69.3 (65.9) 

NA 

ORFlb 

97.8 (98.6) 

NA 

NA 

NA 

82.7 (87.1) 

NA 

ns2 

95.1 (95.0) 

NA 

NA 

NA 

59.0 (59.7) 

60.3 (51.5) 

HE 

96.7 (95.3) 

NA 

89.8 (88.7) 

73.7 (72.9) 

60.8 (58.1) fc 

64.2 (59.5) 

S 

93.5 (91.4) 

92.8 (89.9) 

81.7 (81.0) 

79.3 (79.2) 

66.5 (62.9) 

67.0 (64.4) 

nsl2.9 

96.1 (93.6) 

NA 

89.7 (86.2) 

88.8 (81.7) 

55.0 (47.4) 

59.8 (51.4) 

E 

98.0 (96.4) 

NA 

99.6 (98.8) 

96.1 (91.7) 

72.5 (65.9) 

69.2 (66.7) 

M 

94.8 (94.3) 

NA 

92.1 (93.0) 

91.2 (88.7) 

77.9 (83.5) 

77.5 (82.3) 

N 

96.8 (96.4) 

NA 

94.3 (94.0) 

NA 

71.8 (69.6) 

72.0 (69.9) 

Ia/Ib c 

97.1 (NlT) 

NA 

95.2 (ND) 

NA 

72.3 (ND) 

71.4 (ND) 


a NA, not applicable; the corresponding sequence is not available in the GenBank database. 
b The 5' end of the MHV HE ORF is missing due to a frameshift mutation or sequencing error in NC_001846. 

c The OC43 (ATCC VR759 strain) internal ORF (I) coding region contains a stop codon at position 29345, resulting in two potential coding regions of 60 aa (la) 
and 115 aa (lb). This stop codon is not present in BCoV, which has the capacity to code for a 207-aa protein. This stop codon is also absent in PHEV, MHV, and SADV. 
The percentage of nucleotide similarity is calculated for the continuous Ia/Ib region. 
d ND, not done. 


inoculated into mice (12). The HE gene, another group 2-spe¬ 
cific gene, consists of 1,275 nucleotides and encodes a protein 
of 424 aa. The S gene is located immediately downstream of 
the HE gene and has 4,086 nucleotides. The S protein, which 
consists of 1,361 aa, plays an important role in the attachment 
of the virus to cell surface receptors and induces the fusion of 
the viral and cellular membranes (5, 55). In the genomic region 
between the S gene and the membrane glycoprotein (M) gene, 
two ORFs can be identified: the nsl2.9 gene encodes a putative 
nonstructural protein of 12.9 kDa and is 330 nucleotides long 
(43), while the E gene, 255 nucleotides long, codes for the E 
protein of approximately 9.5 kDa (43). At the 3' end of the 
HCoV-OC43 genome, four major ORFs are present. The M 
gene is 693 nucleotides long and encodes a polypeptide of 230 
aa with a molecular size of approximately 26 kDa (24). The 
membrane glycoprotein is anchored in the viral membrane 
with only a short amino-terminal domain exposed to the exte¬ 
rior of the viral envelope. The nucleocapsid protein (N) gene, 
consisting of 1,347 nucleotides, lies at the 3' end of the genome 
and encodes a 448-aa protein, which is associated with the 
RNA genome to form the nucleocapsid inside the viral enve¬ 
lope. In the 5' part of the HCoV-OC43 N region, two small 
internal (I) ORFs can be identified (la and lb). In BCoV, this 
region is uninterrupted and contains a single I gene which has 
the capacity to code for a 207-aa protein (37). 

HCoV-OC43 sequence similarity to other group 2 coronavi- 
ruses. The sequence similarity among HCoV-OC43, BCoV, 
CRCoV, PHEV, ECoV, MHV, and SDAV was investigated by 
pairwise alignments of the corresponding ORFs and their pro¬ 
teins (Table 2). HCoV-OC43 showed the highest percentage of 
similarity to BCoV in all ORFs except for the HCoV-OC43 E 
gene, which showed 99.6% identity on the nucleotide level and 
98.8% identity on the protein level to the PHEV E gene. 
Maizel-Lenk dot matrix plots illustrate the similarity between 
HCoV-OC43 and BCoV (Fig. 3). 

Phylogenetic analysis. A neighbor-joining phylogenetic tree 
of HCoV-OC43 and 11 other coronaviruses was constructed 
based on an alignment of ORFlb replicase amino acid se¬ 
quences (Fig. 4). As an outgroup, we used the equine torovirus 


(EToV; accession number X52374), belonging to the genus 
Torovirus in the family Coronaviridae. Three phylogenetic clus¬ 
ters, corresponding to the three coronavirus groups, can be 
demonstrated. The SARS-CoV forms a separate branch, al¬ 
though there is strong support for monophyly of SARS-CoV 
with the group 2 coronaviruses, such as HCoV-OC43 and 
BCoV. 

Based on the nucleotide sequence coding for the spike pro¬ 
tein, a maximum-likelihood phylogenetic tree was constructed 
for HCoV-OC43 and several BCoV strains for which the date 
of isolation was known (Table 3; Fig. 5). HEC4408, a corona- 
virus isolated in 1988 from a child with acute diarrhea, was also 
included in the analysis and has actually been shown to be a 
BCoV (66). The time to the most recent common ancestor 
(TMRCA) of HCoV-OC43 and BCoV was dated by three 
methods (Fig. 6). Linear regression of root-to-tip divergence 
versus sampling time situates the TMRCA of HCoV-OC43 
and BCoV in 1891. The maximum-likelihood estimate for 
TMRCA is 1873, with a 95% confidence interval of 1815 to 
1899. The Bayesian coalescent approach dates TMRCA 
around 1890 (95% highest posterior density interval, 1859 to 
1912). This estimate was highly consistent under different de¬ 
mographic models, including an exponential-growth model, 
which resulted in a TMRCA around 1893 (95% confidence 
interval, 1866 to 1918). The evolutionary rate of BCoV was 
also calculated by these three methods (Table 4). A maximum- 
likelihood evolutionary rate of 4.3 X 10~ 4 substitutions per site 
per year was estimated (95% confidence interval, 2.7 X 1CP 4 to 
6.0 X 1CU 4 ). A likelihood ratio test indicated that the molec¬ 
ular clock hypothesis could not be rejected (P = 0.10). 

DISCUSSION 

We report in this paper the first complete genome sequence 
(30,738 nucleotides) of the prototype HCoV-OC43 strain 
(VR759). Until now, only partial sequence fragments of the 
structural protein genes of HCoV-OC43 were available in 
GenBank, leaving the greater 5' part of the genome to be 
determined. The recent discovery of a new human coronavirus, 
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FIG. 4. Phylogenetic analysis of the coronavirus ORFlb replicase 
amino acid sequences. The HCoV-OC43 ORFlb protein (GenBank 
accession number AY391777) was compared to other coronaviruses 
and to an equine torovirus as an outgroup. Group 1, HCoV-229E 
(accession number AF304460), HCoV-NL63 (AY567487), PEDV 
strain CV777 (AF353511), and TGEV strain Purdue (AJ271965). 
Group 2, BCoV strain Mebus (U00735), MHV type 2 (MHV-2; 
AF201929), MHV strain Penn 97-1 (AF208066), and MHV-A59 
(X51939). Group 3, IBV strain Beaudette (M95169), IBV strain LX4 
(AY338732), IBV strain BJ (AY319651). SARS-CoV strain Frank- 
furt-1 (AY291315) is not classified in any of these groups but is most 
closely related to group 2 coronaviruses. Outgroup, equine Berne 
torovirus (EToV; X52374). Regions that were poorly conserved in the 
manually edited multiple protein sequence alignment were deleted 
from the alignment. All columns containing gaps were removed. The 
resulting alignment included 2,083 characters (1,122 being parsimony 
informative) and contained the meld of the following HCoV-OC43 
fragments: 13686-13721, 13737-13793, 13797-13820, 13857-13889, 
13869-13994, 14013-14090, 14127-14174, 14247-14390, 14397-14594, 
14598-14756, 14766-14855, 14859-15230, 15243-15443, 15480-15674, 
15684-15719, 15729-15764, 15786-15854, 15864-15989, 16023-16358, 
16374-16715, 16719-16898, 16902-17093, 17115-17258, 17268-17336, 
17340-17363, 17379-17501, 17535-17561, 17568-17825, 17925-18008, 
18018-18032, 18069-18101, 18207-18440, 18450-18527, 18531-18563, 
18576-18602,18612-18929,18942-19010,19,026-19139,19143-19259, 
19284-19466, 19479-19625, 19686-19793, 20289-20318, 20370-20603, 
20625-20708, 20718-20762, 20769-20885, 20907-21008, 21045-21125, 
21135-21233, 21252-21296, 21309-21431, and 21438-21476. The fre¬ 
quencies of occurrence of particular bifurcations (percentage of 10,000 
bootstrap replicate calculations) are indicated at the nodes. 



FIG. 5. Maximum-likelihood phylogenetic tree of spike gene nu¬ 
cleotide sequences of HCoV-OC43 and several BCoV strains for 
which the date of isolation was known. 


SARS-CoV, necessitates a better understanding of the 
genomic structure and evolution of other known coronaviruses 
in order to gain insights in how this new type could have 
emerged. 

The prototype HCoV-OC43 strain (ATCC VR759) is a lab¬ 
oratory strain that, since its isolation in 1967, has been pas¬ 
saged 7 times in human embryonic tracheal organ culture, 
followed by 15 passages in suckling mouse brain cells and an 
unknown number of passages in human rectal tumor HRT-18 
cells and/or Vero cells. During the passage history, it is likely 
that a number of mutations have accumulated. It would be 
interesting to analyze the complete nucleotide sequence of 
contemporary HCoV-OC43 strains that are free from in vitro 
expansion mutations. 

Nucleotide and amino acid similarity percentages were de¬ 
termined for the major HCoV-OC43 ORFs and those of other 
group 2 coronaviruses (BCoV, CRCoV, PHEV, ECoV, MHV, 
and SDAV). For all HCoV-OC43 ORFs, the highest similarity 
demonstrated was that to the corresponding BCoV ORFs, 


TABLE 3. Date and area of isolation of bovine and human coronaviruses used to calculate TMRCA 


Strain 

Isolation 

date 

Isolation area 

GenBank 
accession no. 

Reference 

BCoV-LY138 

1965 

Utah 

AF058942 

64 

BCoV Mebus 

1972 

Quebec, Canada 

U00735 

18 

BCoV Quebec 

1972 

Quebec, Canada 

AF220295 

32 

BCoV-BECS 

1979 

France 

D00731 

64 

BCoV-COBAAA 

1989 

Giessen, Germany 

M80844 

66 

BCoV-BCQ7373 

1992 

Quebec, Canada 

AF239306 

18 

BCoV-BCQ1523 

1994 

Quebec, Canada 

AF239307 

18 

BCoV-OK05143 

1996 

Kansas 

AF058944 

18 

BCoV-LSU94 

1994 

Louisiana 

AF058943 

9 

BCoV-BC044175 

1997 

Ontario, Canada 

AF239309 

18 

BCoV-OntarioBC043277 

1997 

Ontario, Canada 

AH010241 

18 

BCoV-BCQ3994 

1998 

Quebec, Canada 

AF339836 

18 

BCoV-ENT 

1998 

Texas 

AF391541 

10, 57 

BCoV-LUN 

1998 

Texas 

AF391542 

57 

HEC4408 

1988 

Laubach-Wetterfeld, Germany 

L07748 

William Herbst, personal 
communication 

HCoV-OC43 (VR759) 

1967 

Salisbury, United Kingdom 

AY391777 

42 


Downloaded from http://jvi.asm.org/on March 10, 2015 by guest 





































VIJGEN ET AL. 


J. Virol. 


1602 



0.035 


0.03 3 


0.025 ? 


0.02 


0.015 


0.01 


0.005 


FIG. 6. Results of the evolutionary rate analysis. Line a, linear 
regression of root-to-tip divergence (y axis) versus sampling time (x 
axis). The point at which the regression line crosses the time axis 
indicates the TMRCA (1891). Line b, maximum-likelihood estimate 
(1873) with 95% confidence intervals (1815 to 1899) for the TMRCA. 
Curve c, marginal posterior probability (right y axis) for the TMRCA 
obtained by using the Bayesian coalescent approach. The vertical bars 
in the distribution represent the 95% highest posterior density interval. 
Dates of isolation of FICoV-OC43 and BCoV strains are indicated by 
grey dots. 


except for the HCoV-OC43 E gene, which showed 99.6% iden¬ 
tity on the nucleotide level and 98.8% identity on the amino 
acid level with the PHEV E gene. Based on the high similarity 
between HCoV-OC43 and PHEV in E, and between HCoV- 
OC43 and BCoV in all the other major ORFs, some hypoth¬ 
eses concerning the origin of HCoV-OC43 can be put forward. 
Adaptation of BCoV to a human host and a recombination 
event between BCoV and PHEV leading to a new type of 
coronavirus with a different species specificity could both have 
been responsible for the emergence of a new human corona- 
virus. 

Phylogenetic analysis of coronavirus ORFlb replicase pro¬ 
tein sequences confirms the presence of three coronavirus 
group clusters and a separate branch for SARS-CoV, which 
seems to be most closely related to group 2 coronaviruses (21, 
54). HCoV-OC43 and BCoV cluster together, demonstrating 
the close relationship between the two viruses. There is in fact 
more divergence between the different MHV strains or be¬ 
tween the different IBV strains than between HCoV-OC43 and 
BCoV. The close relationship between HCoV-OC43 and 
BCoV on the genetic level has also been shown to correspond 
to a close antigenic relationship: by using monoclonal antibod¬ 
ies directed against the BCoV S protein, common antigenic 


TABLE 4. Evolutionary rate estimations of the BCoV-HCoV- 
OC43 pair 


Method 

Substitutions/site/yr 

(to- 4 ) 

95% 

Confidence 
interval (10~ 4 ) 

Linear regression 

5.0 


Maximum likelihood 

4.3 

2.7-6.0 

Bayesian inference 

1st codon position 

4.7 

3.1-6.5 

2nd codon position 

3.6 

2.2-4.8 

3rd codon position 

7.8 

5.4-10.4 


determinants for BCoV, HCoV-OC43, and PHEV have been 
demonstrated (62, 63). A phylogenetic tree was also con¬ 
structed for the spike gene of HCoV-OC43 and several BCoV 
isolates for which the date of isolation could be traced. Differ¬ 
ent molecular clock calculations situate the most recent com¬ 
mon ancestor of HCoV-OC43 and the different BCoV isolates 
around 1890. We suggest that around 1890, BCoV might have 
jumped the species barrier and became able to infect humans, 
resulting in the emergence of a new type of human coronavirus 
(HCoV-OC43), a scenario similar to the origin of the SARS 
outbreak. Indisputable evidence for the bovine-to-human di¬ 
rection of the interspecies transmission event, instead of a 
human-to-bovine direction, is not available. However, we con¬ 
sider the occurrence of a 290-nucleotide deletion (correspond¬ 
ing to the absence of BCoV ns4.9 and ns4.8) in HCoV-OC43 
relative to the BCoV genome to be a potential supporting 
argument, as this additional sequence fragment in BCoV is 
also present in MHV and SDAV. Consequently, we assume 
that a deletion from BCoV to HCoV-OC43 rather than an 
insertion in the opposite direction took place during evolution, 
and thus, we hypothesize that the interspecies transmission 
event occurred from bovines to humans. 

Nevertheless, it is possible that two other group 2 coronavi¬ 
ruses, CRCoV and PHEV, might have played a role in the 
emergence of HCoV-OC43. CRCoV appears to be very closely 
related to BCoV and HCoV-OC43 (16), and for the HCoV- 
OC43 E gene, the highest percentage of similarity was found 
with the PHEV E gene, suggesting a possible recombination 
event. To elucidate the evolutionary relationship of HCoV- 
OC43 and BCoV with CRCoV and PHEV, complete genome 
sequence data of CRCoV and PHEV would be required. Mo¬ 
lecular dating has frequently been used to investigate the ori¬ 
gin of viral epidemics (31, 40, 48). The reliability of such an 
analysis is dependent on the validity of the molecular clock 
hypothesis, which assumes that the evolutionary rate is roughly 
constant in the lineages of a phylogenetic tree. Although this 
assumption is frequently violated for viral sequence data (28), 
a molecular clock test indicated that this hypothesis could not 
be rejected for the coronavirus data set investigated here. 

In the second half of the nineteenth century, a highly infec¬ 
tious respiratory disease with a high mortality rate affected 
cattle herds around the world (11, 41). The same disease, or a 
similar disease, is now known as contagious bovine pleuro¬ 
pneumonia (CBPP) and is caused by Mycoplasma mycoides 
mycoides. In the nineteenth century, the clinical symptoms of 
CBPP would have been difficult to distinguish from those of 
BCoV pneumonia, and it can be hypothesized that the bovine 
respiratory disease in the second half of the nineteenth century 
might have been similar to the coronavirus-associated shipping 
fever disease (56). Most industrialized countries mounted mas¬ 
sive culling operations in the period between 1870 and 1890 
(11) and were able to eradicate the disease by the beginning of 
the twentieth century. During the slaughtering of CBPP-af- 
fected herds, there was ample opportunity for the culling per¬ 
sonnel to come into contact with bovine respiratory secretions. 
These respiratory secretions could have contained BCoV, ei¬ 
ther as the causal agent or as a coinfecting agent. 

Interestingly, around the period in which the BCoV inter¬ 
species transmission would probably have taken place, a hu¬ 
man epidemic ascribed to influenza was spreading around the 
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world. The 1889-1890 pandemic probably originated in Cen¬ 
tral Asia (3) and was characterized by malaise, fever, and 
pronounced central nervous system symptoms (53). A signifi¬ 
cant increase in case fatality with increasing age was observed. 
Absolute evidence that an influenza virus was the causative 
agent of this epidemic was never obtained, due to the lack of 
tissue samples from that period. However, postepidemic anal¬ 
ysis in 1957 of the influenza antibody pattern in sera of people 
who were 50 to 100 years old indicated that H2N2 influenza 
antibodies might have originated from the 1889-1890 pan¬ 
demic (45). However, it is tempting to speculate about an 
alternative hypothesis, that the 1889-1890 pandemic may have 
been the result of interspecies transmission of bovine corona- 
viruses to humans, resulting in the subsequent emergence of 
HCoV-OC43. The dating of the most recent common ancestor 
of BCoV and HCoV-OC43 to around 1890 is one argument. 
Another argument is the fact that central nervous system 
symptoms were more pronounced during the 1889-1890 epi¬ 
demic than in other influenza outbreaks. It has been shown 
that HCoV-OC43 has neurotropism and can be neuroinvasive 

(4). 

Maximum-likelihood phylogenetic analysis of the spike gene 
of HCoV-OC43 and several BCoV strains for which the date of 
isolation is known indicates that these strains evolved accord¬ 
ing to a molecular clock. An evolutionary rate on the order of 
4 X 10 -4 nucleotide change per site per year was estimated, 
and this rate was highly consistent across the different methods 
used. This rate falls within the range reported for other RNA 
viruses, including SARS-CoV (14, 50, 51). 

This study provides evidence for viral promiscuity, a phe¬ 
nomenon that has already been reported for several animal 
coronaviruses, including BCoV, for which the potential to in¬ 
fect other species, including humans, has already been de¬ 
scribed (26, 66). The isolation of the SARS-CoV from masked 
palm civets and raccoon dogs indicates that this new type of 
coronavirus was also enzootic in an animal species before sud¬ 
denly emerging as a virulent virus for humans. The character¬ 
ization of the BCoV-HCoV-OC43 pair presented in this study 
provides insights into the process of adaptation of a nonhuman 
coronavirus to a human host, which is important for under¬ 
standing the interspecies transmission events that led to the 
origin of the SARS outbreak. 
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